Indirect Hamiltonian Identification through a small gateway 
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Identifying the nature of interactions in a quantum system is essential in understanding any piiysical phe- 
nomena. Acquiring information on the Hamiltonian can be a tough challenge in many-body systems because 
it generally requires access to all parts of the system. We show that if the coupling topology is known, the 
Hamiltonian identification is indeed possible indirectly even though only a small gateway to the system is used. 
Surprisingly, even a degenerate Hamiltonian can be estimated by applying an extra field to the gateway. 
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I. INTRODUCTION 

When studying any quantum mechanical system, precise 
knowledge of its nature is crucially important. In quantum 
mechanics, any observable phenomena can be explained rig- 
orously, in principle, if we have complete knowledge of the 
system. More specifically, we need to identify the states of 
the system, and the Hamiltonian that governs their dynam- 
ics. Thus, the acquisition of all the relevant information on 
the states and Hamiltonian is essential in understanding how 
nature behaves. The system of interest may include literally 
everything quantum mechanical, from high Tc superconduc- 
tors to microscopic structures in nanotechnology or even some 
highly complex processes in microbiology. 

The full information acquisition is, however, in general very 
hard from an operational as well as from a computational and 
mathematical point of view, even for small systems HI 0. 0] . 
For large many-body systems spectroscopy reveals only little 
information about the Hamiltonian, and generally local ad- 
dressing of its components is required in order to obtain de- 
tails about the system. Spins which can be controlled indi- 
vidually operate as a gateway, through which we can access 
and manipulate the system. A common dilemma is that such 
a gateway not only allows us to interact with the system, but 
also introduces noise to it. From a Hamiltonian identification 
perspective, it is therefore crucial to find minimal gateways 
that suffice to obtain/wW knowledge on the system. While this 
is impossible to answer for generic systems, bounds can be 
derived if the topology of the system is known. In this con- 
text, some positive results have been presented for the case of 
1 -dimensional (ID) chains of spin- 1/2 particles flBIl- That 
is, the coupling strengths between neighboring spins can be 
estimated by accessing only the spin at the end of the chain. 
Since schemes to initialize the state of spins as | || ... J,) by 
operating on the chain end are known |6], such identification 
of the Hamiltonian is sufficient to determine the dynamics of 
the system completely. These results are of interest in their 
own right, yet they were hmited to the simplest of networks, 
i.e., ID chains. 

In this paper, we suggest an estimation scheme for general 
graphs of spins. As well as the details of the Hamiltonian 
identification procedure, we give a precise condition for the 
"gateway" (accessible region) that suffices to make the iden- 
tification possible. For the important cases of finite 2D/3D 
lattices such a gateway is given by one edge or one face of 
the lattice, respectively. This is remarkable because the ra- 




Figure 1: All coupling strengths (black lines) and local magnetic 
fields (blue background) of a 2-dimensional network G — {V, E) 
of spins (white circles) can be estimated indirectly by quantum state 
tomography on a gateway C (enclosed by the dashed red line). The 
coupling strengths and field intensities are represented by the width 
of lines and the depth of the background color, respectively. The la- 
beled spins p and v are used as examples in the proof of the main 
theorem. 



tio between the gateway size and the unknown parameters is 
much higher than in the ID case. We will also show that while 
in the ID case the decay properties of the state in the gate- 
way can identify the Hamiltonian, in the 2D case we need 
its decay properties as well as the transport properties within 
the gateway. Interestingly, our general condition turns out to 
coincide with the criterion for the controllability of spin net- 
works [7]. Our results here thus indicate that Hamiltonian- 
identifiable systems are quantum-controllable and vice versa. 
Furthermore, they support the physical relevance of the topo- 
logical properties discussed later. 

We will study a network with Heisenberg-type interaction. 
This allows us to describe an estimation procedure that is nu- 
merically stable, mathematically simple, and efficient (given 
that we consider arbitrary and large systems). What we at- 
tempt to estimate are the coupling strengths between inter- 
acting spins and the strengths of local magnetic fields. Such 
inhomogeneous fields are very common in experiments, and 
can cause much trouble through dephasing. Hence it is worth- 
while estimating them (such analysis was lacking in JJ, IH). 
Another interesting new aspect we introduce in this paper is 
how to lift degeneracies on the system by applying extra fields 
on the gateway. We show that this is always possible, a result 
which might be relevant beyond the scope of estimation. 
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Our setup is an example of inverse problems that have been 
actively studied in plenty of fields in science and engineering. 
A classical counterpart among those problems that is closest 
to our quantum setting may be the estimation of spring con- 
stants in ID harmonic oscillator chains [8J. However, the res- 
olution to this (classical) problem for generic graphs, even the 
2D case, is still open. It would be intriguing if our results 
in a purely quantum setting could provide some clues to the 
analogous problem in classical settings. 

II. SETUP AND MAIN RESULT 

Suppose that we have a network of spin- 1/2 particles, such 
as the one in Fig [T| We assume that we have knowledge of 
the graph G = (V, E), which describes the network: nodes V 
of the graph correspond to spins and edges E connect spins 
that are interacting with each other. The pairwise interaction 
between spins is Heisenberg type with a known anisotropy A, 
and there is an inhomogeneous magnetic field applied on the 
spins. Then, the Hamiltonian we consider has the form 

(m,n)eE neV 

where Cmn represent the unknown coupling strengths between 
spins m and n, and &„ the unknown intensity of the mag- 
netic field at n, respectively. Here, we also assume Cm„ < 
for all m and n, i.e., ferromagnetic interactions, though the 
setup is readily generalized to other cases. In the above, 
CT*^ {i — X, y, z) are the standard Pauli matrices. The purpose 
of the following will be to estimate c„i„ and fe„ over the entire 
set V of spins by only accessing a small gateway, described by 
a subset C C V (See Fig. [T]i. For almost all practical cases of 
the Hamiltonian identification problem, analyzing the dynam- 
ics in the single excitation sector Hi turns out to be sufficient. 
We will thus denote a single excitation state as |n) S Hi when 
the spin n E F is in the state | f) and all others are in | J,) for 
clarity. The state with all spins in | j) will be written as |0). 

Naturally, the nice challenge here is to obtain information 
about the inaccessible spins C = V\C, which could be the 
large majority of the set. The question is however how small 
can the controlled C be such that we can (in principle) still 
learn all the couplings and fields in V? Intuitively the knowl- 
edge of the graph structure can be useful for making the es- 
timation efficient. For instance, the smaller the number of 
non-vanishing couplings \E\, the more efficiently we can esti- 
mate them. However the efficiency should also depend on the 
structural property of the graph. 

To answer this question, we need to introduce a property, 
known as infecting, of a subset C G V of the nodes JtI M> 
In many-body quantum mechanics this property has 
many interesting consequences on the controllability and on 
relaxation properties of the system (ItISI- The infection pro- 
cess can be described as follows. Suppose that a subset C 
of nodes of the graph is "infected" with some property. This 
property then spreads, infecting other nodes, by the follow- 
ing rule: an infected node infects a "healthy" (non-infected) 



neighbor if and only if it is its unique healthy neighbor. If 
eventually all nodes are infected, the initial set C is called in- 
fecting. The graph in Fig. [T]is an example in which C infects 
V (we encourage the reader to confirm this by coloring the 
nodes in region C and applying the above propagation rule — 
this will make the following proof much more intuitive). With 
this definition, we can summarize the main result of the paper 
as the following 

Theorem: Assume that that C infects V. Then all c„m and bn 
can be obtained by acting on C only. 

This theorem provides an upper bound on the smallest number 
of spins we need to access in order to perform Hamiltonian 
tomography, i.e. given by the cardinality |C| of the smallest 
set C that infects V. To prove the above statement, we first 
present a lemma and its proof. 

Lemma: Assume that C infects V and that all eigenvalues 
Ej (j — I, . . . , \V\) in Hi are known. Assume that for 
all orthonormal eigenstates \Ej) in Hi the coefficients 
(n|£'j) are known for all n E C. Then the c„m and bn 
are known. 

While the assumptions of the lemma may sound unrealistic, 
we will show later how they can be obtained by simple to- 
mography experiments on C. 

Proof of the Lemma: 

We observe that the coupling strengths between spins within 
C are easily obtained because of the relation 

Cnn = (m|i?|n) =Y,Ek{ni\Ek){Ek\n), (1) 

where we defined c,„,„ = (m|iJ|m) for the diagonal terms. 
Since C infects V there is a /i G C and ai/ E C = V\C such 
that 1/ is the only neighbor of fi outside of C, i.e. 

{n\H\fj,) = Vn e C\{iy}. (2) 

For an example see Fig. [T] Using the eigenequation, we obtain 
for all j 

E,\E.f)^H\E,)^Y.^^\^'^H\^)+ E 

mGC neV\C 

Multiplying with and using Eq. (|2]l we obtain 

E,{^l\E,) - cpm(m|i?,) = c^,{u\E,). (3) 

By assumption and by Eq. ([T]i, the left-hand side (LHS) is 
known for all j. This means that up to an unknown constant 
Cpjy < the expansion of \v) in the basis \Ej) is known. 
Through normalization of \v) we then obtain c^^ and hence 
{i'\Ej). Redefining C C U {/i}, it follows by induction 
that all Cmn are known. Finally, we have 

Cmm = (m|i?|m) ^ En~ A ^ c,„„ + 26,„, (4) 
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where N{m) stands for the (directly connected) neighborhood 

of m, and 

Eq = ^ Cmn - X! ^" ^^"^ 

{m.n)£V neV 

is the energy of the ground state |0). Summing Eq. ^ over all 
m and using Eq. (|5]l, we can have the value of X^ney 
thus that of Eq as well, since all other parameters are already 
known. Then we obtain the strength of each local magnetic 
field, bm, fromEq. (|4]i. ■ 



III. TOMOGRAPHY 

Let us now describe how to obtain the information that is 
assumed to be known in the lemma. That is, we need to know 
the energy eigenvalues Ej in Tii and the coefficients (n|£'j) 
for all n g C by controlling/measuring the spins in C. Let 
us first consider the case where the eigenvalues in Tii are 
non-degenerate. The general case will be described in Sec- 
tion UY] To start the estimation, we initialize the system as 
-^(|0) + |1)). As discussed in [6] this can be done efficiently 
by acting on region C only. Then, we perform quantum state 
tomography on the spin n G C after a time lapse t. The entire 
state at t is now 

1 1 / 1^1 \ 

^C/W(|0) + |1)) = ^ le-^^°*|0) + ^/„i(t)|n)J , 

where /„i = {n\U{t)\l) are the elements of the time evolu- 
tion operator in the single excitation subspace. By repeating 
the preparation and tomographic measurements on spin n for 
various times t, we obtain the following matrix elements of 
the time evolution operator as a function of t : 

e'^"*(n|C/(t)|l) =J2{'^\Ej){Ej\l)e-'^^^-^°^K (6) 
j 

If we take n — 1 and Fourier transform Eq. ^ we can get in- 
formation on the energy spectrum of the Hamiltonian in Hi . 
Up to an unknown constant Eq, which will turn out to be ir- 
relevant later, we learn the values of those Ej corresponding 
to eigenstates that have non-zero overlap with We also 
obtain the values of |(l|£'j)p for all eigenstates. Due to the 
freedom in determining the overall phase of a state, we can 
assume that the coefficients for |1) of all \Ej) are real and 
positive, {l\Ej) > 0. Hence observing the decay/revival of 
an excitation at n = 1 we can already learn some Ej and 
all the {l\Ej). This is analogous to the ID case, where this 
knowledge would suffice to obtain the full Hamiltonian [4J . 

In arbitrary graphs however this is no longer the case. In 
fact even if we observed the decay/revival at each n e C we 
would only obtain the | (n|i?j) | ^ , but could not determine their 
phase freely anymore. To obtain the required knowledge for 
the Lemma, we need to observe the transport within C. This 
is represented by Fourier transforming Eq. ^ for n 1, 



allowing us to extract the coefficient {n\Ej) correctly, includ- 
ing their relative phase with respect to {l\Ej). We also obtain 
those eigenvalues Ej which have non-zero overlap with | n) . 
Continuing this analysis over all elements of C, we learn all 
eigenvalues which have overlap with some n G C. Could there 
be eigenstates in Hi which have no overlap with any n £ CI 
The answer is no, as it is shown in |9]. Therefore we can 
conclude that all eigenvalues in the Hi can be obtained. 

Although tomography cannot determine the extra phase 
shift £^0 7 it does not affect the estimation procedure. There 
are three equations that seem to require the explicit values of 
Ej, namely Eq. ([T]i for Cmn inside C, Eq. (O for Cmri and 
coefficients {'n.\Ej) for a spin outside C, and Eq. (|4|i for the 
magnetic fields. It is straightforward to see that for m ^ n 
substituting Ej — Eq into Ej in Eq. ([T]l gives the correct Cmn- 
Similarly, the invariance of Eq. (|4]l is clear as it only depends 
on Ej — Eq. Less obvious is Eq. (|3]l, however, the key is 
that the summation over m G C contains the diagonal term 
c^f, = {fJ'\H\fJ') = £'jI(mI^j>P- Then, by substituting 
Ej — Eq into Ej in the LHS of Eq. (|3]l, it is straightforward 
to confirm that Eq cancels out. Therefore, the precise value of 
Eq is not necessary for the Hamiltonian identification. Even- 
tually, Eq can be calculated by Eq. (|5]l after having all Cmn 
and bn- 



IV. EFFICIENCY AND DEGENERACY 

The efficiency analysis of the Hamiltonian tomography is 
roughly the same as in |i4J. Due to the conservation of exci- 
tations, the sampling can be restricted to an effective di- 
mensional Hilbert space, and the speed is some polynomial in 
|V^|, provided localization is negligible. One difference how- 
ever is that in arbitrary graphs it might be less likely that the 
spectrum is non-degenerate. An explicit example can be given 
for a square lattice with equal coupling strengths, with the 
spectrum Ekj — Ek + Ej, k,j — 1,...,N, where the Ek 
the ID energies of the corresponding chain. A uniform ID 
system on the other hand would typically be non-degenerate. 
Of course "exact degeneracy" is highly unlikely; however ap- 
proximate degeneracy could make the scheme less efficient. 
Here, we suggest to lift degeneracies by applying extra fields 
on the gateway C. Since C is only a small subset of the spin, 
it is not obvious at all that this is possible. We prove the fol- 
lowing perhaps startling property of the infection property: 

Theorem: Assume that C infects V. Then there exists an op- 
erator Be on C that lifts all degeneracies of H in the 
single excitation subspace. 

Proof: 

We will prove the above by explicitly constructing a Be that 
does the job. This Be will be very inefficient and even re- 
quires full knowledge of the Hamiltonian, but is only intro- 
duced here for the sake of this proof. Let us denote the 
eigenvalues of H as Ek and the eigenstates as where 
d = 1, . . . , D{k) is a label for the I? (A;) -fold degenerate 
states. Let us first concentrate on one specific eigenspace 



4 



I l-E^), d = 1, . . . , D(fc)} corresponding to an eigenvalue Ek- 
Since the eigenstates considered here are in the single excita- 
tion subspace, we can always decompose them as 

\Et)cc = \4>i)c ® \0)c + |0)c ® \^t)c, (7) 

where we introduced the unnormalized states \4>'i)c and 
\'4>t)c in the single excitation subspace of C and C, re- 
spectively. As shown in [9] we know that \4>'l)c 7^ Vd. 
This is because if there was an eigenstate in the form of 
|0)c ® \'4>t)c '^hen applying H repeatedly on it will necessar- 
ily introduce an excitation to the region C, in contradiction to 
being an eigenstate. In fact the set { l^^)^, d — 1, . . . , Z3(fc)} 
must be linearly independent: for, if it was linearly de- 
pendent, there would be complex numbers akd such that 
'}2d^kd\<P't)c — 0, and because the eigenstates are degen- 
erate, i2d^kd\E'l,)cc = J2dO^kd\0)c \tPt)c would be 
an eigenstate with no excitation in C, again contradicting 
Ref. H- This leads to an interesting observation that the de- 
generacy of each eigenspace can be maximally |C|— fold, be- 
cause there can be only |C| linearly independent vectors at 
most in the single excitation sector on C. Thus minimal in- 
fecting set of a graph gives us some bounds on possible de- 
generacies. 

Now we consider a Hermitian perturbation B^c ® Ic ('^^ 
be specified later) on the system and compute the shift in ener- 
gies. We shall see that it suffices to assume that Bfec 1 0)c — 0. 
In first order, we need to compute the eigenvalues of the per- 
turbation matrix 

cc{Et\Bkc\Et')cc =c {cf^i\Buc\cf^i)c. (8) 

Can we find a Bkc such that all eigenvalues differ? For 
that, note that {\4>f)c^ = 1, ■ ■ ■ , D(k)^ are linearly inde- 
pendent, which means that there is a similarity transform Sk 
(not necessarily unitary, but invertible) such that the vectors 
l^^)^ = 5*^7^10^)0 are orthonormal. The perturbation ma- 
tiix can then be written as c{(,t\^kEkcSk\£,t )c- If we set 

SlBucSk^Y.^kd\ek)c{ii\ 

d 

we can see that the Hermitian operator Bkc = 

i:d^kd(sl)^\ek)cmS^' gives us energy shifts 
ekd- Therefore, as long as we choose the €kd mutually 
different from each other, the degeneracy in this eigenspace 
is lifted by Bkc- This happens for an arbitrarily small 
perturbation A^. We choose Afe such that the lifting is large, 
but in a way such that no new degeneracies are created, 
i.e. IIAfeSfccll < min (Ai?) , where AE* are the energy 
differences of H. However, the perturbation XkBkc may 
well lift other degeneracies of H "by mistake". Note that 
by construction Bkc conserves the number of excitations in 
the system (See Eq. ([8]l). Therefore, we can now consider 



the perturbed Hamiltonian H' = H + XkBkc and find its 
remaining degenerate eigenspaces in Hi. Naturally, the 
number of degeneracies with H' is less than that with H. 
Following the above procedure, we pick one eigenspace 
and find an operator Bk'c that lifts its degeneracy. Keeping 
IIAfc'Sfe'cll < min(Ai?') we continue to add perturbations, 
until we end up with a sum of perturbations that lift all 
degeneracies in TYi. ■ 

The above theorem demonstrates that degeneracies can in 
principle be lifted. In practice, we expect that almost all op- 
erators will lift the degeneracy, with a good candidate being 
an inhomogeneous magnetic field on C. One could even ran- 
domly choose operators on the gateway C until the system 
shows no degeneracies. Albeit being inefficient, our theorem 
shows that this strategy will eventually succeed. Note also in 
the theorem it sufficed to consider operators within Hi, i.e. 
E = Em.nec^™"l™)("l with b*^^ = 6„„, so maximally 
|Cp parameters need to be tested. For instance, if the system 
is a chain, B necessarily corresponds to a magnetic field on 
spin 1. 



V. CONCLUSIONS 

We have shown how a small gateway can efficiently be used 
to estimate a many-body Heisenberg Hamiltonian, given that 
the topology of the system is known. It is surprising to see 
how a simple topological property of a network of coupled 
spins - infection - implies so many far-reaching properties, 
from control to relaxation, from the structure of eigenstates to 
possible degeneracies, and, as we have shown here, for Hamil- 
tonian identification. 

Our results can be seen as an example of inverse problems 
in quantum setting. It would be intriguing to explore a possi- 
ble link between ours and similar problems in classical setting, 
such as 2D graphs of masses connected with springs. Also, it 
would be interesting to study if the methods of [5], which does 
not require state preparation, can be applied to this setup. A 
further application could be found, for example, in estimating 
the hidden dynamics in an environment of an controllable sys- 
tem, such as a nanoscale device |12]. Of course, generalizing 
the present results to a wider class of many -body Hamiltonian 
will be important from both theoretical and practical perspec- 
tives. 
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